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ABSTRACT 

' Are magnetic fields important in primordial star formation? Assuming that star formation occurs via 

I an accretion disk that is turbulent, initially because of local gravitational instability, we calculate the 

, disk structure for realistic accretion rates. We predict that local gravitational viscosity is able to drive 

^ ' accretion, without the disk fragmenting. We then estimate the rate of dynamo amplification of seed 

O ! magnetic field. Turbulence in a stratified disk can be helical, with different signs of the helicity in each 

hemisphere. This provides a key ingredient for production of global scale magnetic fields whose sign of 
flux is sustained over many orbit times. The resulting fields can drive coUimated protostellar outflows 
that reduce the star formation efficiency from the initial gas cloud, especially once the protostar has 
contracted to the main sequence, at ~ lOOM©. We estimate that the outflows are powerful enough to 
eject some material from the host dark matter halo and to initiate relatively strong magnetization of 
^ . the local intergalactic medium. Close to the protostar, the outflow acts to shield the disk and equatorial 

\jr^ ' regions from radiative feedback, such as ionizing photons, and this may enable accretion up to relatively 

If^ , large stellar masses. We conclude that magnetic fields cannot be ignored from models of primordial star 

■ formation. 

Subject headings: early universe — MHD — stars: formation — stars: magnetic fields — stars: winds, 
i outflows 

m 

O ' 1. INTRODUCTION 

Primordial star formation has attracted much attention recently, due to advances in numerical models of cosmic structure 
formation that can resolve down to solar system scales (Abel, Bryan, & Norman 2002, hereafter ABN), constraints on the 
Q ■ reionization epoch from the Wilkinson Microwave Anisotropy Probe observations of the microwave background (Kogut 
I et al. 2003), and the discovery of very metal poor Galactic halo stars (Aoki et al. 2002; Christlieb et al. 2002), whose 
c/j . elemental abundances may be the flngerprints of the earliest supernovae. 

' The typical mass of the first stars is an important, but unknown, quantity. It may be set by feedback processes local to 
the forming star. To investigate this possibility it is necessary to model the star formation process by including as much 
physics as possible, whilst making allowance for a realistic geometry of the accreting gas. Magnetic fields play crucial roles 
in present-day star formation, but have usually been assumed to be unimportant in the primordial case. It is the goal of 
this paper, which builds on the detailed, yet non-magnetic, star formation model of Tan & McKee (2003a), to critically 
examine this assumption. 

At typical locations in the standard cosmological model, baryonic over-densities collapse to form the first stars at 
redshifts ~ 20 (ABN, Bromm, Coppi, & Larson 2002). These first "complex" baryonic systems, consisting of a protostar 
and accretion disk, live for many dynamical times so that amplification of seed magnetic field by dynamo action may be 
possible. If the field strength becomes dynamically significant, then an outfiow can be generated. 

The ubiquity and energy of molecular outflows from present-day protostars (Zuckerman, Kuiper, & Kuiper 1976; Kwan 
& Scoville 1976; Edwards & Snell 1982; Bally & Lada 1983) was a surprise and forced a radical change in theoretical models 
of star formation. Outflows are now thought to play a key role redistributing angular momentum in the accretion disk 
(Blandford & Payne 1982), setting the star formation efficiency (Levreault 1984; Myers et al. 1988; Matzner & McKee 
2000), and depositing energy into the star- forming interstellar medium (Norman & Silk 1980). It is often speculated 
that magnetic fields and outfiows may be unimportant for primordial star formation because the strength of seed fields 
is presumably very weak. Thus they have generally been ignored in studies of the collapsing gas cloud (e.g. Haiman, 
Thoul, & Loeb 1996; Tegmark et al. 1997; Nakamura & Umemura 1999; ABN; Bromm et al. 2002) and in models of 
the formation of the protostellar core (Stabler et al. 1986; Omukai & Nishi 1998; Ripamonti et al. 2002; Tan & McKee 
2003a). However, as in the local Universe, the presence or absence of dynamically-important magnetic fields and outflows 
from primordial protostellar disks has important implications for the basic star formation process. 

The ampliflcation of seed magnetic field by dynamo action in Population III protostars and its ejection in outfiows may 
also bear on the long-standing question of the origin of galactic and intergalactic field strengths (e.g. Kulsrud et al. 1997, 
and references therein), particularly if this mode of star formation was widely and evenly spread throughout the Universe. 
Such a model would be analagous to those which propose that fields are ejected from active galactic nuclei (AGN) (Daly 
& Loeb 1990; Colgate & Li 1999; Kronberg et al. 2001), though the AGN models would produce more ordered flux on 
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large scales. If field amplification in Population III star formation is very efficient, then much of the intergalactic magnetic 
field may have been in place before or at the same time as the era of reionization, which has implications for models of 
field generation in cosmic ionization fronts (Subramanian, Narasimha, & Chitre 1994; Gnedin, Ferrara, & Zweibel 2000). 

Starting from conditions representative of the end state of the numerical simulation of ABN, Tan & McKee (2003a, 
hereafter TM) predicted the mass infall rate, accretion disk structure, and protostellar evolution associated with primordial 
star formation. Using these results, in this paper we consider a wider range of possible accretion disk models, calculate 
their structure over the evolution of the protostar and assess their gravitational stability with respect to fragmentation 
(§2). We discuss the ability of these disks to generate large scale magnetic fields via a dynamo (§3), which is likely to be 
a key requirement (Lubow et al. 1994; Blackman 2003; Blackman & Tan 2003) for magnetically mediated outfiows from 
thin disks (e.g. Blandford & Payne 1982: Lovelace et al. 1987; PcUetier & Pudritz 1992; K5nigl & Pudritz 2000). We 
show that strong, dynamically important fields can be produced even from the weak seed fields expected of primordial 
protostellar disks. We address the implications of feedback from hydromagnetic winds on the star formation process in 
§4, and conclude in §5. 

2. PRIMORDIAL PROTOSTELLAR ACCRETION DISKS 

The dark matter halos and gas clouds relevant to primordial star formation form from hierarchical mergers of smaller 
structures, including the accretion of matter along filaments (ABN). Thus in general we expect that the initial conditions 
for star formation are gas cores with non-zero angular momentum. The collapse should proceed to a disk, which then 
channels mass to the protostar via viscous torques. 

We use TM's model of the rotating, freely falling accretion envelope that develops inside the sonic point, rgp, of the 
collapsing gas. The rotation is parameterized by 

„ ^ t'rot(?-sp) ^ ^rot(rsp) 

■^"""^ ~ t^Kep(r.p) (GM,p/r,p)i/2' ^> 
which is the ratio of the mass- weighted rotational velocities, Vrot, relative to Keplerian, WKop, of infalling matter at a radial 
distance from the protostar corresponding to the mean location of the sonic point. Inside this region, which contains a 
total mass Mgp, we assume conservation of angular momentum (Ulrich 1976) so that matter falls freely to the midplane 
and forms a disk of size 

= /Lp^P ^3.44 (1^)' (i±fi)"' (I|)"V'-»AAU. P) 

For the numerical evaluation we have used the mass-radius relation expected at the sonic point of the collapsing cloud 
(TM), which involves the entropy parameter, K' , of the polytropic equation of state of the cloud, relative to the fiducial 
case with T = 300 K at a hydrogen density of 10^ cm""^ — larger values of K' correspond to denser gas cores. We have also 
assumed an efficiency, c^^d, of collapse from the core to the disk, which is unity when feedback processes are unimportant, 
and related Mgp to the stellar mass, m*, at the time when Mgp has collapsed to the disk: m*(l + fd) = Mgp, where fd is 
the disk mass relative to the star's and with the assumption that there is little diversion of the accretion flow from Vsp to 
Td- From the results of the simulations of ABN we have taken a fiducial value of /kop = 0.5, but there is likely to be a 
significant dispersion depending on the formation history of the halo and cloud. The fiducial value of fd adopted by TM 
is 1/3, because the disks are expected to be relatively massive since their angular momentum transport must initially be 
driven by gravitational instabilities. 

Material falls onto the disk (and directly to the star) at all radii r < Vd (Figure 2 of TM; Cassen & Moosman 1981). 
The characteristic disk scale is much greater than the stellar radius (r* < lOOi?0 = 0.5 AU, TM) so most matter accretes 
to the star via the disk. The evolution of r* as the stellar mass increases is shown in Figure la, together with the linear 
growth of Td- Although infall onto the disk provides some effective viscosity, this is not important in the very inner regions 
and so we assume a standard "agg" model for the viscosity (Shakura & Sunyaev 1973). We discuss the choices for the 
value of ttss below. 

The overall rate of collapse and thus accretion rate of the star is given approximately by m* ~ m^/tg, where tg is 
the local free-fall timescale evaluated at the density of the region that contains mass m* at the moment it undergoes 
dynamical collapse. The accretion rate is thus set by the structure of the pre-stellar gas core, which in turn is set by 
the cooling properties of trace amounts of molecular hydrogen in the primordial gas. A more careful evaluation of the 
accretion rate to the star (TM) yields 

= 0.026(1 + fd)-''/'eTK'^'/' yr-\ (3) 

These high accretion rates make it likely that the disk will build itself up to a mass significant compared to the stellar mass 
(we have set fd = 1/3) and will become susceptible to gravitational instabilities. Two dimensional simulations of clumpy, 
self-gravitating disks show self-regulation with a^s — {iltth)~^ up to a maximum value ttgs — 0.3 (Gammie 2001), where 
n is the orbital angular velocity, <th = SfcT'c.d/(o'7off d) the thermal timescale, S is the surface density, Tc.d is the disk's 
central (midplane) temperature, and Tog-.d the effective photospheric temperature at the disk's surface. Fragmentation 
occurs when ntth ^ 3: this condition has the best chance of being satisfied in the outermost parts of the disk that are 
still optically thick (see below). In addition to local gravitational viscosity, global (m = 1 mode) instabilities (Adams et 
al. 1989; Shu et al. 1990) may be efficient at driving inffow if the disk becomes too massive. 
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In addition to gravitational instabilities, the magneto- rotational instability (MRI) may become important, yielding 
viscous stresses that correspond to somewhat smaller values of ass- This process could be particularly important in the 
inner accretion disk, where gravitational instability is suppressed. Values often quoted for the MRI are Ugs ^ 0.01 (Balbus 
& Hawley 1998), but larger values have also been measured and neither simulations nor physical understanding have yet 
determined why it cannot differ by an order of magnitude (or whether the constant formalism is really appropriate 
for MRI-driven turbulence). We evaluate the ability of the MRI to operate in §3. If the disk is stable with respect to the 
MRI, then ass may be set by molecular viscosity and thus much smaller. In this case the accretion rate may effectively 
go to zero. However, the disk is still being fed at a rate given by equation (3), so that the surface mass density inevitably 
builds up until self-gravity (and the associated larger values of ttss) becomes important. 

Given these considerations, we calculate the radial structure of the inner disk at any given point in the protostcUar 
evolution by assuming it is fed smoothly at a rate given by equation (3). We also use the standard theory of steady, 
thin, viscous accretion disks, with a spatially constant viscosity parameter, ass, and ignore energy injection from the star. 
The viscosity is assumed to be a function of the total pressure, though it is gas pressure that dominates over radiation 
pressure in the earlier stages (m* < 20Mq), which we think are relevant for field growth (§3). We evaluate cases with 
ass = 0.01,0.3. 

The inner scale of the disk is set by r^, , which is predicted from the protostellar evolution model of TM. It is somewhat 
sensitive to the choice of ass for m* < 20Mq because this sets the temperature of gas in the inner disk and this thermal 
energy is advectcd into the star (Fig. la). We adopt primordial gas abundances X = 0.76 and Y = 0.24 so that in the 
atomic phase /x = 1.22toh and nne = 0.079nH, and use the opacities of zero metalhcity gas of Rogers & Iglesias (1996) 
for T > 6000 K and of Lenzuni, Chernoff, & Salpeter (1991) for T < 6000 K (however, here we are interested primarily in 
the higher temperature regime). The method assumes that the disk's central temperature, Tc^d is much greater than the 
surface temperature, Tcff^d, i-e. that the disk is very optically thick: this condition is not well satisfied during the very 
earliest stages, or for the outer disk, since the opacity is set by and declines rapidly for T < 7000 K. 

We calculate the ionization state of H and He from the Saha equation and include the energy needed for ionization of 
H and He and the thermal energy in the disk's energy equation: 
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where F is the flux emitted from one side of the disk, and eth and e/ are the thermal energy and dissociation/ionization 
energy per unit mass, averaged over the thickness of the disk. The first term on the right hand side of eq. (4) is the 
viscous dissipation per unit disk face area. Note that usually den^/dr < and dlj/dr < 0. The efi'cct of the thermal 
energy term is small for thin disks such as those we consider, so we drop it from our numerical calculations. However, the 
dissociation/ionization term, which is often neglected, can be an order of magnitude larger than thermal term and can 
have a major impact on the disk structure. We therefore retain this term. We follow Frank, King, & Rainc (1995), using 
an approximate solution based on radiative diffusion so that F ~ AaT^^^/Sr, where Ted is the midplane temperature and 
T is the optical depth from the midplane to the surface. 

Dividing the disk into many discrete radial zones, we solve the full set of disk structure equations (e.g. Frank et al. 
1995), starting from the outer regions where the gas is almost completely neutral. The high accretion rates of primordial 
protostars lead to the ionization energy being quite important, particularly in the carlicT stages of evolution before the 
protostar starts contracting towards the main sequence. For complete ionization of atomic gas, the energy absorption rate 
is 

dE\ o ( rh* \ , , 

^ m,ei,„a. = 2.4 x 10=^ L©, (5) 



calculated assuming a dissociation energy of 13.6 eV for H ionization, and 24.6 eV and 54.4 eV for the first and second 
stages of He ionization.^ 

Figure 2 shows three examples of disk structure for ags = 0.01, with to* = 1, 10, IOOMq. At these masses the stellar sizes 
are taken to be r* = 100, 300, 4 i?© and the accretion rates (eq. 3) are m* = (17,6.4,2.4) x 1O~^M0 yr~^, respectively. 
Figure 3 shows the equivalent models with a^s = 0.3. 

The self-gravity in the disks can be gauged by considering the Toomre stability parameter for a Keplerian disk, Q = 
Cg^l {'kGYi), which is < 1 in unstable regions. The outer parts of the disk are most susceptible to gravitational instability. 
Following a line of reasoning similar to that presented by Goodman (2003) concerning self-gravity in quasar accretion 
disks, consider a disk with ass = 0.01 provided by the MRI: the model will be physically self-consistent inside the radius at 
which Q = 1. At larger radii, we expect that gravitational viscosity will become important, allowing larger values of ass so 
that Q remains approximately equal to imity. The results of Gammie (2001) show that once agg ^ 0.3 then fragmentation 
occurs. This corresponds to the location in our a^g = 0.3 models where Q = 1. From Figure 3 wc can sec that Q > \ 
for all the regions of the disks with Tc,d ^ 7000 K, and so fragmentation should not occur in these regions. At larger 
radii we expect the disks to become optically thin because of the temperature dependence of H~ opacity (approximately 
cx T^^"'), and have thermal timescales long compared to the orbital period. This appears to occur before Q < 1, so that 
fragmentation would then not occur anywhere in the disk. 

^ The results of ABN indicate that the first solar mass or so of material to be accreted would be almost fully molecular, which for a dissociation 
energy of 4.48 eV for H2 increases the coeflacient of equation (5) by 11%. We do not include this minor correction factor in our numerical 
calculations. 
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Fig. 1. — (a) Upper panel: Evolution of the radius of the protostar, r*, and the radial extent of the ionized (fi{H) > 0.5) region of the 
accretion disk, rjon for models with viscosity parameter Ogs = 0.01,0.3 (dashed, dotted lines, respectively). The formation is from the collapse 
of a core with K' = 1 and /xep = 0.5 (see TM for more details). The full extent of the accretion disk (r^) is shown by the long-dashed line. 
(b) Lower panel: Mean number of orbits of fluid elements in the ionized region of the disk for Ogs = 0.01, 0.3 (dashed, dotted lines). The gap 
at m* ~ IOM0 is when rion = ft because of the expansion of the protostar during this stage of the evolution (TM). 



The infall of material onto the disk provides an effeetive viscosity and source of turbulence that is likely to be quite 
important in the outer disk, and that this will also help prevent fragmentation. Even if Q became < 1, fragmentation 
is not inevitable: global modes of gravitational instability, such as spiral waves, could develop instead. These modes, if 
present, are relatively efficient at transporting angular momentum and could then connect the outer infall region with 
the inner, viscous accretion disk, where first local gravitational viscosity, and then perhaps viscosity from the MRI, bring 
material all the way to the immediate vicinity of the protostar. The final accretion to the stellar surface may be mediated 
in part by angular momentum loss in protostellar outflows, the properties of which are discussed in §4. 

We conclude that fragmentation of typical primordial protostellar accretion disks does not occur, if the results of 
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Gammie's (2001) simulations are applicable. However, further numerical studies are desirable to confirm this prediction. 
We note that in the most advanced simulations to date (ABN) of the global primordial star formation process, which 
probe conditions in the initial collapse down to scales of order astronomical units (just before formation of a hydrostatic 
protostellar core and an optically thick accretion disk), there is also no indication of fragmentation. Primordial star 
formation appears to favor single (or perhaps binary) stars rather that star clusters. The reservoir of gas available to 
the forming star is several hundreds to thousands of solar masses, but the typical mass of the stars is likely to be set by 
feedback from the protostar acting on its own accretion. Following TM, here we shall consider the uninterrupted growth 
of the star at rate given by eq. (3), accreting from a stable accretion disk, and use this as the framework for assessing 
magnetic field generation and subsequent feedback effects. 

3. DYNAMOS AND PROTOSTELLAR DISKS 

Stabler et al. (1986) and Omukai & Palla (2001, 2003) considered the structure of primordial protostars in the limit 
of spherical accretion. They found that the protostar is radiatively stable (apart from a very thin outer layer, ^ 10~^ of 
the stellar mass; Omukai 2003, private comm.) until m* ~ lOOM©, when nuclear H burning creates a convective core. 
Thus, except at very high masses, little dynamo amplification in the protostar is expected. We therefore study helical disk 
dynamos, focusing on the production of large scale ordered fields because these are the most susceptible to escape into 
the corona (Blackman 2003; Blackman & Tan 2003). There they can drive outflows that produce a feedback on the star 
formation process. A non- helical dynamo may also be present, but this produces disordered random fields that do not 
drive a strong outflow. For the helical dynamo to operate, turbulence in a stratified rotator is the crucial requirement. 

3.1. Basic Theory 

We discuss likely sources of turbulence, the ability of MRI-driven turbulence to produce large scale flelds, and then the 
helical dynamo. 

3.1.1. Sources of Turbulence 

We expect that the initial driver of turbulence is gravitational instability (Gammie 2001). This may then be followed 
by MRI-driven turbulence (Balbus & Hawley 1998). 

For sufficiently ionized and magnetized disks, the combination of radially decreasing angular velocity in a rotating 
plasma and seed magnetic field can in principle drive the MRI (Balbus & Hawley 1991, 1998). The vertical stratification 
then provides the turbulence with the helical property from which a large scale helical dynamo feeds. However, the critical 
wave number below which the MRI grows, depends on the initial field strength: there is a minimum initial field for the 
instability to drive turbulence in the disk. Since primordial star-forming environments have only very weak seed fields 
(discussed later), it is important to evaluate this critical field strength. It is determined by requiring that the MRI growth 
time is shorter than the diffusion time at the scale of the initial critical vertical growth wavenumber associated with the 
initial seed field (Balbus & Hawley 1991;1998). This wavenumber is /c^ = 3Q/va- Since the growth rate is of order f2, the 
condition that a given dissipation process is slower than the growth rate requires 

xk^jn = 3xf)M < 1, (6) 
where x is a particular microphysical diffusion coefficient to be considered. Note that if x is an ambipolar diffusion 
coefficient, then no constraint on the field emerges since then x v\. However, if x is the Spitzer magnetic diffusivity 
(Spitzer 1967), then the condition becomes 

R\2 / X 1/2 / \ -3/2 /, .^ 1/2 

-D \ . _ 4^-3/2 I rn^ \ ' I r \ '/InA^' 



where A is the usual Coulomb logarithm and T3 = T/IO^K. 

From this estimate we can see that while the MRI process may be able to operate given the expected initial conditions 
of present-day star formation, in the early universe the fields are likely to be too weak: we estimate in §3.2 that the 
initial fields that thread the disk have a maximum value of order 10~^^ G. Yet, even if the MRI does not work initially, 
gravitational instability in the disk can in principle provide local turbulence that is amenable to an ags treatment (e.g. 
Gammie 2001; Wada et al. 2002). As we argued in §2, lacking other sources of viscosity, the accretion disk should 
establish a state in which the accretion is driven by gravitational instability. The angular momentum transport could 
also be global, via modes that are not describable with the ags formalism. The issue of which mechanism operates and 
at what levels is not fully resolved (Balbus & Papaloizou 1999), requiring further numerical study. Here we assume that 
local turbulence does arise. 

In this context and for the interpretation of simulations, it should be noted that the meaning of ass is sometimes 
misunderstood: a turbulent disk incurs large fluctuations so that simulation data must be averaged temporally or spatially 
for appropriate comparisons to the ags formalism to be made. The formalism for the helical dynamo that we discuss below, 
only requires that the stratified helical turbulence yields a viscosity that can be modeled by an ass parameter in the mean 
statistical sense. 

— 1/2 

In a turbulent disk, the small scale field grows faster than the helical large scale field by a factor of order ass 
(Blackman 1998, 2003). Thus, well before the large scale field growth saturates, the small scale field will grow to the 
critical value required for the MRI to operate. For this reason, even if the initial turbulence is driven by gravitational 
instability, we expect the later phase of field amplification to be mediated by the MRI and its associated ass- 
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Fig. 2. — Protostellcir disk structure for models with ass = 0.01 and m, = 1,1O,1OOM0, for which r* = 100, 300, 4ilQ and rh* = 
(17,6.4,2.4) X 10~^Mq yr~^, respectively. From top to bottom the panels show (1) surface density, S; (2) ratio of scaleheight to radius, 
hjr, and ratio of gas pressure to total pressure, /3; (3) midplane ionization fractions of H+, He+, He'^+; (4) disk midplane temperature, ^ 
(the dotted lines show results for when the ionization energy is neglected); (5) number of orbits, rtorbi Toomre Q stability parameter, and 
Rosseland mean opacity k, evaluated at the midplane. Note that all quantities axe azimuthal and temporal averages of the disk, which, being 
turbulent, exhibits local fluctuations. 



3.1.2. Comment on MRI and Large Scale Fields 

Field amplification by the MRI in a non-stratified disk is sometimes said to produce large-scale fields (Balbus & Hawley 
1998), but the meaning is subtle and easily misunderstood. In the simplest MRI simulations that use a non-stratified 
periodic box, the saturated magnetic energy spectrum emerges as Kolmogorov with a kr^/^ shape. No global flux can 
grow in a periodic box. The toroidal field has the largest coherent scale, but this is only of order the disk height. Most 
importantly, although the energy of this field is sustained, the sign of the field fluctuates rapidly, on a timescale of order 
the turbulent correlation time (about equal to a Keplerian rotation time for MRI unstable disks). The reason is that 
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Fig. 3. — As Figure 2, but with ass = 0.3. 



scales of order the box thickness are themselves unstable to the MRI, and thus are the largest turbulent scale. The largest 
scale fields produced by the non-stratified MRI have the same spatial and temporal properties of the largest scale of 
the turbulence (Hawlcy ct al. 1995). Such turbulent scale fields can be shredded before buoyantly rising to the coronae 
(Blackman 2003; Blackman & Tan 2003), and are therefore not ideally suited for driving outflows. 

However when stratification is present, the MRI can produce helical turbulence. "Helical" describes the property that 
rising turbulent eddies in the northern (southern) hemisphere expand and twist clockwise (counterclockwise) to conserve 
angular momentum, while falling turbulent eddies twist counterclockwise (clockwise). Both rising and falling plasma 
maintain the same sign of pscudoscalars such as (v • V x v) and (b • V x b) , (where v and b are the turbulent velocity and 
turbulent magnetic field in Alfven units, respectively), with each having opposite signs in opposite hemispheres. Such 
finite pseudoscalars enable amplification of global scale fields by helical dynamo action when the rising plasma is threaded 
by a magnetic field (e.g. Moffatt 1978; Parker 1979; Brandenburg 2001; Blackman & Field 2002). In a real non-periodic 
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system the helical dynamo can maintain a global large scale flux over many rotation periods, as seen in non-periodic box 
simulations of disks (Brandenburg et al. 1995; Brandenburg & Donner 1997). Note that Stone et al. (1996) claimed that 
no helicity was observed in their stratified periodic box disk simulations, when averaged over hemispheres (Stone 2003, 
private comm.). Further studies are required to reconcile these results. 

3.1.3. The Helical Dynamo 

Recently, Blackman & Field (2002) developed a dynamical nonlinear helical dynamo theory that accounts for the 
nonlinear backreaction of the magnetic field on large scale field growth and agrees with recent helical dynamo simulations 
(Brandenburg 2001). The driver of exponential large scale field growth, aa, turns out to be proportional to the residual 
helicity (v • Vxv) — (b • Vxb) (Pouquet, Prisch, & Leorat 1976), rather than just the kinetic helicity (v • Vxv) used in 
textbooks for similar dynamos. Large scale field growth occurs in two phases: an early, fast kinematic phase, and a late, 
slow resistive phase. Generation of the large scale field is associated with large scale magnetic helicity, whilst helicity 
conservation means that the small scale helicity grows with the opposite sign. As this happens, (b • Vxb) grows and 
offsets the (v • Vxv), shutting down the dynamo. These principles seem to survive in dynamos with shear (Blackman 
& Brandenburg 2002). Although the above results are for simplified dynamos, we shall use the saturation principles in 
what follows, by incorporating them into an a — 17 accretion disk. Other approaches to mean field accretion disk dynamos 
have been developed by, for example, Pudritz (1981), Riidiger, Elstner, & Stepinski (1995), Stepinski (1995), Reyes-Ruiz 
& Stepinski (1997; 1999) and Campbell & Gaunt (1999). 

To estimate the large scale field, B. produced by a helical dynamo, we need the growth rate. We use a mean field 
theory for large scale field growth in a turbulent disk. This approach presumes that the mean quantities vary on time 
scales long compared to that of the fluctuating turbulent quantities. Therefore, the derived mean field growth rate, 7, 
should be slower than the eddy turnover rate, i.e. 7 < f2, where the latter is also the growth rate of the random field (e.g. 
Balbus & Hawley 1998). We shall find that this condition is indeed satisfied. 

We assume B = 'BoExpl'-ft + ik^z], where is the vertical wavenumber. The basic equations for the radial and toroidal 
field in an "a^^" dynamo for a Kelperian disk become (e.g. Ruzmaikin et al. 1988) 

jBr = -ik, [adB^) - (3klBr (8) 

and 

jB^ = -{3/2)nBr - pklB^, (9) 

where we ignore spatial dependencies of the pscudoscalar a^i and turbulent transport scalar 0, and assume averages arc 
taken within a given hemisphere. We assume ad = qctdo where the time-dependent multiplicative factor q < 1 represents 
the backreaction from (b • Vxb), with q = 1 in the kinematic regime. 

There are two possible dominant contributions to ado- a Goriolis force effect and a magnetic buoyancy effect. First 
consider the Goriolis force acting on the largest turbulent eddies, of typical scale I. A rising and expanding (falling 
and contracting) eddy rotates in the opposite (same) sense as that of the imderlying disk rotation during its turbulent 
correlation time, ted = I /v. Therefore, while conserving angular momentum about the disk axis, such eddies gain a specific 
angular momentum about their own symmetry axis such that (v • Vxv) is non- vanishing and of opposite sign in each 
hemisphere. Using the equation for Vxv evolution and the mass continuity equation, and assuming that the density 
varies on scale h, this leads to the estimate (Ruzmaikin et al. 1988) 

ado = tcd(v • Vxv) ~ tl^{v ■ (Vx {n X v))) ^ ^(v . flV • v) ~ — x Mm[Rb, 1], (10) 

where the last relation follows from mass continuity, and Rb = is the Rossby number. 
For thin disks in the Shakura-Sunyaev formalism, 

lr.al/'h/Rl/' (11) 

and 

v^all^Rl^^, (12) 

which follow since 1/ vl v'^/{Rbfl) = a^sCsh and O/i ~ Cg. Using (11) in (10), we have 

ad = qadn ^ qa^^^h. (13) 
It also follows that the maximum fractional kinetic helicity satisfies 

h = Min[l, ado/v ^ all^/Rl'\ (14) 

It is important to note that while equation (10) is a standard estimate for adn, Brandenburg & Donner (1997) found that 
for large scale field generation in local shearing box simulations, the contribution seemed to have the opposite sign of that 
required to produce the measured mean field. It was later suggested (Brandenburg 1998) that a magnetic buoyancy term 
of the form ted{t>xby{gz/p)) where bx,by are the fluctuating components of the magnetic field, is the vertical gravity, 
and p is the pressure, actually has the correct sign. Since this term is ^ {l/^){v\/cl)ass{v'^/r){h/r) ^ assCs ^ flP/h, 
which has the same order of magnitude and scalings as equation (10), it does not significantly affect our calculations based 
on the use of equation (13), unless the two cancel to high order. The cycle period analyzed by Brandenburg & Donner 
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(1997) suggests some cancellation. However, they studied the late time evolution of the dynamo, after the kinematic 

regime saturates. 

With our assumed form for B, we solve equations (8) and (9), use equation (13), and set kz ~ 1/h and (3^1/ = assCgh, 
to derive the growth rate 

7 = (3arfA;,fi/4)V2 _ ^^2 ^ c.V'c«(x/3g/2 - aU^)/h. (15) 

Note that with q = 1, for growth this requires a dynamo number D = \{3/2)flaQh^ / > 2 = Dc, the critical dynamo 
number. Critical dynamo numbers for disks are sometimes derived to be > 8 (e.g. Ruzmaikin et al. 1988) but this value 
depends on the assumed profile of a{z). We have crudely assumed a nearly constant profile with ^ > 0, and thus our 
critical dynamo number is low. 

At a given stage in the evolution when the stellar mass is m* , we estimate B as the smaller of (i) the saturation value 
at the end of the kinematic regime and (ii) the value reached given the growth time that the disk dynamo has been able 
to operate, 

where tinit is the time at which the protostellar disk becomes susceptible to turbulence. We define f = at the time when 
the protostellar core has just started to form, then from equation (3) we have 

t(m,) = 27(1 + /,)i°/^e;;°/^;^'-i5/^(m,/Mo)i°/7 yr. (17) 

Several factors potentially determine Unit- Firstly, we require a critical level of ionization so that the field's growth time 
due to the dynamo is substantially shorter than the diffusion time. Brandenburg et al. (1995) found a critical fractional 
ionization of disk material of ^ 10^^** for growth of the MRI, corresponding to a neutral-ion collision frequency ^ 10 
times the orbital frequency. The MRI has a growth time ass times that of the large scale helical field and so our critical 
ionization fraction is ~ a~^10~^'^ ~ 10"^'^. This condition is well-satisfied for the inner regions of primordial protostellar 
disks, once the star has a mass that is even just a small fraction of a solar mass. Thus tinit (and the corresponding 
TO^jnit) is set by the time for disk formation in the collapsing flow, i.e. when the centrifugal radius of the flow becomes 
significantly larger than the protostellar radius. For the flducial case of K' = 1 and fuep = 0.5, this occurs very early on, 
when the star has mass m*,init — OAMq. 

One might think that the time available for fleld growth at a given radius is the accretion timescale at that radius, 
rather than the disk age. However, since the mean field formalism is statistical, amplification occiirs as long as the system 
remains in quasi-statistical equilibrium (i.e. as long as helically turbulent plasma is steadily supplied to that loc;ation). A 
related point is that as long as (3 > [L^ where L ^ h in the vertical variation scale of the large scale field, a poloidal 
mean field is likely not frozen into the plasma on an accretion timescale and would diffuse away if the helicity effect 
is ignored. This then motivates the need for the aj, effect and large scale dynamos in the first place (Lubow et al. 1994; 
Blackman 2003; Blackman & Tan 2003). For completeness however, we will also provide the estimated field when the 
growth time is limited to the local accretion time. 

After the kinematic regime, subsequent field growth becomes resistively limited as the kinetic and current helicities 
nearly balance. The magnitude of the helical component of the mean field when this occurs is given by (Blackman & 
Field 2002; Blackman 2003) 

Bh ~ (47rpA)V2^,(z//j)V2 ^ (47rp)V2a^^c„ (18) 

where the latter similarity follows from equations (11), (12), and (14). We take this to be an estimate of the poloidal 
fleld, i.e. we assume Br ~ Bh. We assume that the strength of the toroidal field at the surface B^, is comparable to 
its value within the disk. This toroidal fleld is larger than Br because it is linearly amplifled above the value of Br hy 
shear as it rises into the coronae during a buoyancy or diffusion timescale. These timescales are thuo ~ /i(47ryo)^/^/B^ 
and idif = h/assCs, respectively. Linear growth of B^f, above the level predicted by equation (18) in a time tbuo gives 
B(p = Br^t\yuo = ass''^Cs(47rp)^/^, which from equations (18) and (12) implies Br/B^ = ali^ . Similarly, if the growth of 
B^ is limited by turbulent diffusion, then B^ = Brfltdif, so that Br/B^ ^ agg. 
To summarize: for an initial seed fleld Bq we expect the flnal field strengths to be 

~ Min{BoExp[7oigrow]> assCs(47rp)^/2} , (19) 

where 70 — aW^Cg/h = aW^Q, is the kinematic 7, obtained from equation (15) with 5 ~ 1, and 

B^ = a-' Br (20) 

if B^/(47r/9)^/^ < UssCs, (i.e. when thuo > idif) or 

B^ = a-y^Br (21) 
if B^/{4TTpy/^ > ttssCs, (i.e. when thuo < idif)- The latter case applies in the saturated state. 
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3.2. Expected Values for Star-Forming Disks 

To determine the expected B in star-forming disks, we must estimate the seed field, assess flux freezing for the initial 
collapse of the cloud to form the disk, and estimate the growth timescale to determine if the field has saturated. The 
seed field for the primordial star-forming environment can be estimated from the Biermann battery mechanism for non- 
barotropic flows. On_scales much larger than that of an indvidual star, numerical simulations indicate that by ^; = 18, 
seed fields of order Bq ^ 10~^^ Gauss can be generated (Kulsrud et al. 1997). We assume such a field is generated 
inside the pre-stellar core region in the simulation of ABN, which has a size ^ 100 pc and a density ~ 10^^® g cm^"^. The 
infall speeds within this region are of order kni/s, vastly greater than the drift speeds of the ions through the neutrals 
(the ionization fractions are much greater than the equilibrium values for the gas, which has a minimum temperature of 
^ 200 K). Thus field is frozen into the contracting gas: going from parsec to AU scales the field strength increases by 
factors of ^ 10-^", so the seed field in the disk may then be ^ 10^^^ G. 

It may actually be that much larger seed fields could arise in the protostar from the Biermann battery term if the pressure 
and density gradients are misaligned inside the protostar as well. The protostellar Biermann battery (e.g. Kulsrud et 
al. 1997) value reached after a rotation time (at which point other growth terms become important) can be estimated 
to be Bo ~ (mHc/e)^^^ ~ 2 x 10-^\T/10^Ky^^{r/10^^cm)-\h/r) Gauss, where we have used fl = cjh and have 
conservatively taken V 1/r. 

In the disk, fiux freezing and field advection are not guaranteed since, as discussed above, turbulent diffusion of the 
mean field may be relatively efiicient and the accretion timescale may be relatively long. Thus, in addition to a fiducial 
initial seed field strength of Bq = 10^^*'G, wc also consider a conservative value of Bq = 10^^*'G. Since the disk-dynamo 
amplifies the initial field exponentially, the uncertainties in the strength of the seed field do not significantly affect the 
time to reach saturation. 

Using 7o ^ ali'^Cs/h = all^il from (15) with q = 1, and approximating tgrow = t{mt,) = 0.7m*/m*, the first term on 
the right of equation (19) becomes 

3/2 , ^ , _o/2 / • \ -1 



(22) 



while the second term gives 

Br ~ 49G(a3s/O.Ol)i/2(m,/Af0)i/4(m,/O.OlM0yr)i/2(r/lOi3cm)-3/4(/i/loi2cm)-V^ (23) 

where mass continuity at r = 2r*, p = rh*/(47r/i^assCs), was used, and we choose parameters commensurate with the 
results of §2. Although primordial protostars are thought to have accretion rates somewhat higher than the 1O~^M0 yr~^ 
adopted in equation (22), it is clear from these equations that wc expect the field to saturate when the protostar has 
a relatively small mass. Note that the seed fields for present-day Galactic protostellar disks would be many orders of 
magnitude higher than the primordial case: the initial Galactic interstellar material possesses a field of > 10~^ G. Thus 
for a wide range of parameters equation (23) is the relevant estimate for the dynamo- amplified large scale field strengths 
in present-day, as well as primordial, protostellar disks. For the primordial case, we evaluate the saturation point more 
accurately in the to* versus plane, and show the results in Figure 4. 

In the fiducial case the field in the inner disk saturates by about IMq, and this does not depend very sensitively on aga- 
in lower angular momentum cores the disk emerges later, at which point the dynamo rapidly causes the field to saturate. 
Although we expect that mean field dynamo can grow field at a given location for longer than the local accretion time, 
we also evaluate the cases where tgrow = ^acc ^ 27rnorb/^ — r'^ / (agsCsh) . The resulting field strength must always be 
smaller than that from equation (22) for disks for which Ugs models are self-consistent. The saturation point is delayed 
until masses m, ~ 10 — 2OM0 in the fiducial case. 

4. IMPLICATIONS OF LARGE SCALE FIELDS FOR PRIMORDIAL STAR FORMATION 

4.1. The Power of Magnetically Mediated Outflows 

The primary effect of the magnetic field that we consider is the generation of a magnetically mediated jet or wind. 
Integrating the vertical Poynting flux through the disk in each vertical half gives an estimate of the magnetic luminosity 
available for driving each outflow: 

Lma,g,l/2 ~ B^^sBz,s^r^ , (24) 

where B^^g and Bz.s are the toroidal and vertical fields just inside the disk surface (defined at one scaleheight above 
the midplanc). The field strengths wc estimated in the previous section represent fields inside the disk, so the use of 
equation (24) really requires a solution of the global boundary value problem (e.g. Riidiger, Elstner & Stepinski 1995; 
Reyes-Ruiz & Stepinski 1997;1999). For present purposes, we make crude estimates for the surface fields. We take 

^ -B^.avo = ctss^^"^ Br.nvc SiS estimated in the previous section, where -B^.avo ^ (3/4)-B0 and Br^a^c ^ (3/4)_Br are 
the respective average values over a scaleheight, using flux freezing in which the density falls a factor of 1/e between the 
midplane and surface. Note that for a constant ass disk, va/cs is a constant (Blackman 1998). Then with Cg oc r~^/^ and 
p cx r~^/^, we have Br oc r~^/'^, so that imag is dominated by the contribution generated at the inner disk radius. 

Now we use the above relations to estimate Bz,s by appealing to V • B = 0. We integrate this over an annular wedge of 
the disk extending from the midplane to the surface, with azimuthal extent A0, and with radial extent 5r, where r 
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Fig. 4. — Conditions for saturation of the disk magnetic field. Protostellar growth starts from low masses and with accretion rates that 
decline with time and mass: the fiducial case is shown by the straight thick solid line. The rising curves represent the condition that the disk 
magnetic field has been able to saturate to values given by equation (23), when the disk is old enough to have allowed sufficient amplification 
of the seed field of 10~^^ G. The disk dynamo is assumed to start operating at i^it = O.'lAfg, which corresponds to when the centrifugal 
radius of the accretion fiow is about an order of magnitude greater than the stellar radius in the fiducial case with /xep = 0.5. The dashed, 
solid, and dotted lines show the saturation condition for ass = 0.001, 0.01, 0.1, all with r = lO'^^cm and m*_init = 0.4Mq. With ass = 0.01, the 
long-dashed line shows this condition for r = 10^* cm, the dashed-long-dashed line shows it for a seed field of 10"'^® G, and the dot-long-dashed 
line shows it evaluated with m*_init = 4.7Mq, corresponding to an initial core with /xep = 0.1. If the dynamo amplification time is limited by 
the accretion timescale then saturation occurs at masses indicated by the short line sections (dashed for ass = 0.01 and dotted for ass = 0.3), 
which use the results for riorb of §2 and assume a seed field of 10~^^ G. Note that these results for n^vh are only valid for m* and rh* close 
to the fiducial case shown by the thick solid line. 



is the characteristic radial scale over which the vertical and radial fields change significantly. The toroidal flux does not 
contribute to this integral, and by the assumption that the vertical field vanishes at the midplane, no net fiux penetrates 
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the equatorial plane. Then integrating V • B = over the described volume gives 
A(j) 



R h - fi t nY^^ fh 



= A<^\ \B.A'\rT - rT)\ , (25) 



(OL \ ^ TTXi 77% f T \ — ^ 

mJo.OIM^ vr-i ITo13^J ^^^^ 



where the subscripts i and o refer to inner and outer radii, and we assume hjr = hi/vi, that is, the ratio is independent of 
r, which is a reasonable approximation (see Figs. 2 and 3). In the limit rg '> Ti and using Br,ave = (3/4)5^ from above, 
we obtain 

B,,ilBr,i = (9/16)(/i/r-). (26) 
Prom equation (21), B^^s = Oi^^^ Br,a.ve = (3/4)a^^^^i?r,i so that equation (24) for the outflow power becomes 

i^mag,i/2 ~ (27/64) (Vr)a-'^'^r^^'- (27) 

The procedure for using equation (27) to assess the value of Lmag for star-forming disks, is to pick a particular star 

formation model (which means choosing a^^, iris, and to*) and inferring from its solution whether the field has saturated 
or not (from eqs. 22 and 23). The appropriate field strength is then used in equation (27) to estimate the wind luminosity. 
The case of a saturated field (eq. 23) is the most interesting: 

ir„ag,i/2 - (27/64)al/2^^ = Q.MaTL^^^y^, (28) 

1/2 

Mq O.OIM0 yr- 

whcrc ij„ax.i/2 — G'TO*TO*/2r is the upper limit on the magnetic luminosity through each disk plane available from 
accretion. Note that imag,i/2 < -^max,i/2 as expected, for all reasonable values of a^s (formally ags < !)• The evolution of 
i-mag with m* is shown in Figure 5a, together with the radiative luminosity from the fiducial protostellar evolution model 
of TM. The wind luminosity is a fixed fraction of the accretion luminosity from the inner disk. 

Fields produced by a helical dynamo arc ordered such that their smallest gradient scale exceeds that which would 
otherwise result in their turbulent destruction before rising to the corona. But there is an additional process in the 
generation of global large scale fields in the corona that we have not discussed above: the fields supplied by the disk 
dynamo must incur further evolution to even larger scales once in the corona. Their scale in the disk is a few times h, 
but the needed large scale in the corona is ^ r. While this further opening or "dynamical magnetic relaxation" process 
requires further study with direct application to the disk-corona system, there is ample reason to believe that it happens 
sufficiently fast when the supplied fields from the dynamo are helical (as they would be from a helical disk dynamo, and 
further supplemented by foot point shear): It has long been known that helical fields relax to the largest scale available 
to them in magnetically dominated environments (Taylor 1986) and it has been further shown recently that when helical 
magnetic fields are injected into such a magnetically dominated region, this relaxation occurs on an Alvfcn crossing time 
(Blackman & Field 2003). The speed is independent of the value of the magnetic diffusivity, as long as the latter is 
not identically zero. For this reason the dynamical relaxation in the corona would naturally allow further opening of 
disk-supplied fields to global scales, not unlike the "magnetic carpet" evolution in the solar corona (Schrijver & Zwaan 
2000). 

Given the magnetic luminosity, we estimate the mass-loading and momentum of the outflow. We assume that at 
large distances from the star, the energy content of the flow is dominated by the kinetic energy of the matter, so that 
Lraa-g = T^^v^l^, whcrc TO^ is the mass outflow rate and is the terminal velocity. Theoretical models of magneto- 
centrifugally driven outflows suggest that is about equal to the escape velocity from the protostar, t^esc,* (e.g. Shu et 
al. 2000; Konigl & Pudritz 2000). We write = fvVesc,*, where /„ is a factor of order unity. Then the mass flux of the 
outflow is 

27 a'^^^ 

rriu, = -^-^'rn* ^ 0.084m*. (30) 

We note that ass is uncertain, as are the numerical prefactors of equations (29) and (30). For a massive protostar that is 
approaching the Eddington limit, the mass-loading of the outflow may be significantly enhanced. 

4.2. Feedback Effects from Outflows 

In present-day star formation, protostellar outfiows contribute the following feedback effects: (1) they prevent the 
accretion of gas in the polar regions of the star-forming core that would otherwise have joined the star and they divert a 
fraction of the accretion flow in the disk (e.g. Matzner & McKee 2000, hereafter MM); (2) they inject mechanical energy 
into giant molecular clouds, which helps to keep the clouds turbulent (Norman & Silk 1980); (3) through the presence 
of dense, outflowing gas near the star, they inhibit the reach of other feedback effects, such as ionizing radiation (Tan & 
McKee 2003b). We discuss the significance of these processes in the context of primordial star formation. 
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4.2.1. Star Formation Efficiency 



First, consider the expulsion of gas from the polar regions of the core. By balancing magnetic tension and magnetic 

pressure gradients, Matzncr & McKcc (1999) (sec also Shu ct al. 1995; Ostrikcr 1997) showed that far from the star the 
angular distribution of the momentum in a radial hydromagnetic wind can be approximated by 



where /i = cos9 and is a small angle, which is estimated to be ^ 0.01 for winds from low-mass stars (Matzncr & McKec 
1999), but may be somewhat larger for winds from massive stars. We adopt this value in the primordial case (our results 
are not particularly sensitive to this choice). Equation (31) is also a reasonable description for the angular distribution of 
the mass flux in the wind, assuming the terminal velocity (far from the star, but before deceleration by the ambient gas) 
is independent of 9. 

To evaluate the efficiency of star formation from a core (i.e. the fraction of mass that forms a star) Ccorei we follow 
MM and find the angle, ^esc = cos~^/iesc, where the wind sweeps up core material to the surface escape speed of the core, 
Wesc.c- The analysis assumes thin, radiative shocks, purely radial motion and monopole gravity. We consider cores with 

angular mass distributions of the form dM/dQ = {1/4:Tt)Q{h)M with Q(/x) = (1 — ^^)"/ (1 — ^^)" dfx, and evaluate 
models with n, ranging from 0, isotropic, to 4, which describes a flattened distribution that mimics the effect of some 
rotational support. In the notation of MM, the escape condition is given by (1 + 9^ — fil^^)Q{^csc) = mf/{XM), where 
X = 0.132cg[ln(2/(?o)/ln200]iJcsc,c,5[(l5ti)/"^*)/40kms~^]~-'^, and Cg is a factor of order imity that accounts for the effects 
of gravity on shock propagation. Using the results of MM, we estimate Cg ~ 4.6, for steady winds that decouple from 
the swept-up shell at the core edge in a core with kp = 20/9 (which is the fiducial value of TM). These cores have 
Wcsc.c = 3.22(M/1000AfQ)^^/^if '■''/^ km/s. Wc ignore the influence of material beyond the core "edge", which is partly 
balanced by our neglect of the wind's influence on the shell beyond this point. The efficiency is given by 



in equation (3), this definition of efficiency includes the effect of diversion of material from the disk into the outfiow. 
However, since /„, = m^/m* < 1, Ccore — e*d- 

In summary, wc arc assuming that the properties of the protostellar outflow, and how they evolve with m*, do not 
depend on ecorc- In reality, as ecorc and m, become smaller, then the protostellar evolution changes, e.g. the stellar growth 
time becomes longer relative to the Kelvin-Helmholz time, so that the star would tend to be smaller. However, since Pw 
is integrated over the stellar evolution, it should not be so sensitive to these changes. 

We show the evolution of pw/m* with m* for our fiducial model in Figure 5b. Initially r* is approximately constant and 
the specific momentum of the wind grows only because of the increasing mass and escape velocity from the star. Then 
above m» ~ 20Mq, the star starts shrinking towards the main sequence and there is a jump in the momentum input. 
This marks the start of the transition to lower efficiency star formation from the core (Fig. 5c) . This transition may have 
important implications for the final mass of the star, and will be examined in more detail in a future paper. Note that 
this efficiency is not particularly sensitive to the parameter 6o when the critical angle for ejection of gas from the core is 
> ^0- That is, mass limits to star formation resulting from this process do not depend on Bq so long as Oq <C 1. 



We have seen that the protostellar outflow is eventually able to break free of the core. With the distribution of 
momentum given by equation (31), this is particularly true near the rotation axis, since about half of the initial outflow 
mass is concentrated inside 6o- A lOOM© star releases a wind with total momentum of about 10^ Mq km/s, with about 
a quarter of this inside 9q of each jet, for which the injected mass is about 2Mq. Each central jet region can then sweep 
up ^ 500A'/q to about 5 km/s, which is the approximate sound speed and escape speed of the gas in the vicinity of the 
dark matter mini-halo. The fiducial core of TM has a total virializcd baryonic mass of ~ lO^M© (based on the simulation 
of Abel et al. 2002), with a few xW^Mq of gas still infalling within a couple of hundred parsecs. For reasonable values 
of 6o < 0.1, and even with allowance for turbulent entrainment of matter at the sides of the jets and lateral spreading, 
the outflows should easily be able to eject some gas from the dark matter halo. Thus we expect the extent of mechanical 
feedback to reach beyond the immediate environment of the star-forming core. 

The physical implications of this feedback are difficult to assess, since radiative fciedback from the massive protostar 
is also acting on the cloud during this time. If outflows do escape from the halo, then they will start to pollute the 
intergalactic medium with their magnetic fields. This is only one possible mechanism for producing such fields. Others 
include field generation in shocks caused by supernovae or ionization fronts (Subramanian et al. 1994; Gnedin ct al. 2000), 
and these mechanisms should have a longer reach than the protostellar outfiows. Outfiows from AGN are also likely to 
be important, although these may be more sparsely distributed than Population III stars. 



dpw Pw 1 



(31) 



dfl 47rln(2/6io)(l + 6i^ 





4.2.2. Feedback on the Larger-Scale Cloud and Break-Out from the Dark Matter Halo 



4.2.3. Confinement of Radiative Feedback by Protostellar Outflows 
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Fig. 5. — Evolution of protostellar luminosities, cumulative specific outflow momentum, and the resulting star formation efficiency for the 
fiducial /kcp = 0.5 and K' = 1 case, with ass = 0.01. (a) Top panel: Evolution of wind luminosity, Z/mag, is shown by the dashcd-long- 
dashcd line. Also shown arc the total bolomctric luminosity from the protostar (solid line) including contributions from the iimcr (r < lOr*) 
accretion disk (dotted line), stellar accretion luminosity including from direct spherical accretion and boundary layer accretion (dashed lino), 
and internal protostellar luminosity (long-dashed line) (see TM for more details). The Eddington luminosity is shown by the dot-long-daslicd 
line, (b) Middle panel: Cumulative momentum of the protostellar outflow divided by the stellar mass. This estimate assumes that the wind 
speed is the escape velocity from the protostellar surface and that at large distances from the star the energy content of the flow is dominated 
by kinetic energy, (c) Bottom panel: Efficiency of star formation due to erosion of the gas core by the outflow. The density profile of the 
initial core is specified by dM/dCl = (l/47r)Q(//)M, with // = cos6l and Q(/i) = (1 - /i^)"/ jj(l - ^^^n ^^j^ gou^ ling 

is n = (isotropic core), 

dotted is n = 1, dashed is n = 2, long-dashed is n = 3, dot-dashed is n = 4. 



The presence of dense, outflowing gas near the protostar can suppress the influence of radiative feedback. For example, 
the propagation of ionizing photons should be impeded. The characteristic density at the base of the flow region near the 
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star is 

- " ^ (t) " (bt) ( 10-3;;,,- ) (it^fc) (loS;) ('^) 

where /„, = rn^/rh^. and wc have assumed the wind is launched from the inner part of the disk from an annulus of inner 
radius /^r* = 2r«, and outer radius 

Assuming an isotropic flux, the ionizing luminosity necessary to balance recombinations and fully ionize the gas, is 
given by 

-^crit f'^^"''' 
^'^ Jfrr, 

where a^^^ = 2.62 x 10~^^ s~^ is the radiative recombination rate to the excited states of hydrogen at 10^ K and rig is the 
electron number density. We assume singly ionized He, so that 

s„,. = 2.9 X 10-/.- {^y (liF^lj^)^ (t^j;)" <^^> 

Models of Population III stellar structure (Tumlinson & Shull 2000; Baraffc, Hcger, & Woosley 2001; Bromm, Kudritzski, 
& Loeb 2001; Schaerer 2002) yield zero age main sequence (ZAMS) ionizing luminosities that can be approximated via 
S = 9.42 X 10''^(m,/30MQ)2-5phs-i for (12 < m^/M© < 40), S = 9.17 x 10'*9(m,/100Mo)i-^phs-i for (40 < tu^/Mq < 
120) and S = 3.75 x lO5°(m*/3OOM0)i-2 ph s"^ for (120 < m^/MQ < 1000). The critical ZAMS mass that can ionize the 
outflow from a protostar accreting at the fiducial rate (Eq. 3) is then given by 

/ f \ -0.23 / f \ 0.46 

m.,eHt = 34/-0-^'^ i^t^ (^^ j i^'0-9«Mo, (36) 

so long as 12 < m*,crit/AfQ < 40. Similar relations can be derived for when m*,crit is in a higher mass range. 

Primordial protostars only rc^acli the ZAMS at masses around IOOM0 (Omukai & Palla 2001; 2003; TM). Before this 
the ionizing luminosity is initially very small because the protostar is large and relatively cool. In the fiducial case the 
ionizing luminosity, including the contribution from accretion luminosity, then increases rapidly to become equal to the 
ZAMS value at about 40Mq (Figure 7 of TM). Therefore, we expect the outflow to be able to confine ionizing radiation 
for much of the pre-main sequence phase of protostellar evolution up to about 4OM0, but not after that. These results 
have implications for feedback models that arc based on the ionization of the protostellar disk (e.g. HoUenbach et al. 
1994). In present-day massive star formation this process of confinement by the outflow can create extremely compact 
H II regions around massive protostars (Tan & McKee 2003b). A more complete analysis of the feedback effects in the 
primordial case, including those from FUV photons, is beyond the scope of this paper and is deferred to a future study. 

5. CONCLUSIONS 

We have presented calculations of the radial structure of primordial protostellar disks, using realistic models for the 
rate of accretion and the size of the protostar (Tan & McKee 2003a). We have followed the evolution of disk properties as 
the protostar grows in mass and changes in size. These disks differ from those around present-day protostars, because of 
the lack of cooling due to dust and because of the high rates at which they are fed matter from the collapsing core. The 
high accretion rates and large protostellar sizes require inclusion of the ionization energy in the structure equations, the 
effect of which is to reduce the disk temperatures and luminosities, particularly for disks around protostars with masses 
< 2OM0. We argue that fragmentation in the disks is unlikely to occur, if local gravitational instabilities can provide an 
effective viscosity of the magnitude seen in the simulations of Gammie (2001). 

Wc have shown that in situ dynamos can produce dynamically significant large scale magnetic fields in primordial star- 
forming disks. For present-day disks, such fields would likely be important for driving outfiows in all star-forming systems. 
For the primordial case, the disks of the lowest mass protostars would not have enough time to generate strong fields, 
but the majority of the star-forming parameter regime at higher masses leads to strong flelds even for very conservative 
assumptions about the initial seed magnetic field. Thus models of primordial star formation should include the dynamical 
infiuence of magnetic fields. 

We emphasize that for these conclusions to be valid, gravitational instability, rather than the MRI, must provide the 
initial source of turbulence in the disk. Even though the ionization fraction is high, the seed field is below the critical 
value required at the mode of fastest growth. This contrasts with the situation in present-day protostellar disks. A second 
requirement is that the turbulence in the disk must generate helicity. It is important for this stage of the evolution to be 
investigated with numerical simulations that can resolve the helicity. 

If large scale magnetic fields can be amplified to their saturation strength, then their main effect is to drive a collimated 
wind. This reduces the star formation efficiency by diverting material from the accretion disk and by sweeping-up ambient 
gas from the polar regions, that otherwise would have collapsed. The star formation efficiency is significantly reduced once 
the protostar reaches about IOOM0. The outfiows are probably strong enough to eject matter from the dark matter mini- 
halo and to start magnetizing the intergalactic medium. The presence of outfiows also increases the density distribution 
of gas near the star in the region above the accretion disk. The full effects of this have not yet been worked out, but it 
appears likely that the propagation of ionizing photons is inhibited, at least while the protostar is in the earlier stages of 
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contraction towards the main sequence. These results will be applied to a more detailed study of protostellar feedback in 
a forthcoming paper. 
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